# Load the 'plotrix' library
library(plotrix, lib = './lib/')

coords = read.csv("scratch/hohmann.coords", header=FALSE)[,c(2,3)]
other.coords = read.csv("visualization-data/hohmann.1004.coords",header=FALSE)[,c(2,3)]

xy = rbind(coords, other.coords)
plot(xy, asp=1, main="Hohmann transfer orbit: scenario 1004", xlab="x-coordinate", ylab="y-coordinate")

# Draw a big green earth
draw.circle(0,0,6357000,nv=1000,border="green",col="green",lty="solid",lwd=1)

# Draw initial orbit radius
draw.circle(0,0,6457000.000001785,nv=1000,border="black",col=NA, lty="solid", lwd=1)

# Draw target orbit radius
draw.circle(0,0,38330909.090909086,nv=1000,border="red",col=NA, lty="solid", lwd=1)
